taDUfcO 3 © DDT] 
ENDOCRINOLOGY 



ORIGINAL RESEARCH ARTICLE 

published: 14 February 2013 
doi: 10.3389/fendo. 2013. 00004 



Association between polymorphisms in the aryl 
hydrocarbon receptor repressor gene and disseminated 
testicular germ cell cancer 

Leon J. S. Brokken 1 *, Yvonne Lundberg-Giwercman 1 , Ewa Rajpert-De Meyts 2 , Jakob Eberhard 3 , OlofStahl 3 , 
Gabriella Cohn-Cedermark 4 , Gedske Daugaard 5 , Stefan Arver 6 and Aleksander Giwercman 1 - 7 

' Department of Molecular Reproductive Medicine, Lund University, Malmo, Sweden 

2 Department of Growth and Reproduction, Rigshospitalet, Copenhagen, Denmark 

3 Department of Oncology, Skane University Hospital, Lund, Sweden 

4 Department of Oncology-Pathology Radiumhemmet, Karolinska Institute and University Hospital, Stockholm, Sweden 

5 Department of Oncology, Rigshospitalet, Copenhagen, Denmark 

6 Centre for Andrology and Sexual Medicine, Karolinska University Hospital Huddinge, Department of Medicine, Stockholm, Sweden 

7 Reproductive Medicine Centre, Skane University Hospital, Malmo, Sweden 

In the Western world, testicular germ cell cancer (TGCC) is the most common malignancy 
of young men. The malignant transformation of germ cells is thought to be caused by 
developmental and hormonal disturbances, probably related to environmental and 
lifestyle factors because of rapidly increasing incidence of TGCC in some countries. 
Additionally, there is a strong genetic component that affects susceptibility. However, 
genetic polymorphisms that have been identified so far only partially explain the risk of 
TGCC. Many of the persistent environmental pollutants act through the aryl hydrocarbon 
receptor (AHR). AHR signaling pathway is known to interfere with reproductive hormone 
signaling, which is supposed to play a role in the pathogenesis and invasive progression 
of TGCC. The aim of the present study was to identify whether AHR-related poly- 
morphisms were associated with risk as well as histological and clinical features of TGCC 
in 367 patients and 537 controls. Haplotype-tagging single-nucleotide polymorphisms 
(SNPs) were genotyped in genes encoding AHR and AHR repressor (AHRR). Binary 
logistic regression was used to calculate the risk of TGCC, non-seminoma versus 
seminoma, and metastasis versus localized disease. Four SNPs in AHRR demonstrated a 
significant allele association with risk to develop metastases (rs2466287: OR = 0.43, 
95% CI 0.21-0.90; rs2672725: OR = 0.49, 95% CI: 0.25-0.94; rs6879758: OR = 0.27, 
95% CI: 0.08-0.92; rs6896163: OR = 0.34, 95% CI: 0.12-0.98). This finding supports the 
hypothesis that compounds acting through AHR may play a role in the invasive 
progression of TGCC, either directly or through modification of reproductive hormone 
action. 
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INTRODUCTION 

Testicular germ cell cancer (TGCC) represents the most-frequently 
diagnosed cancer in young men between 20 and 40 years of 
age. The rates vary considerably according to region and ethnic- 
ity, which suggests that both environmental and genetic factors 
contribute to TGCC. In Europe, a three- to fourfold increased inci- 
dence has been noted during the past 30-40 years, which suggests 
that an increased exposure to persistent environmental pollutants 
may play an important role. 

Common persistent environmental pollutants include dioxins 
and polycyclic aromatic hydrocarbons, which exert their toxic 
and carcinogenic effects through the ligand-activated transcrip- 
tion factor aryl hydrocarbon receptor (AHR). In complex with 
its binding partner AHR nuclear translocator (ARNT), it medi- 
ates cellular responses to xenobiotic compounds. The AHR/ ARNT 



heterodimer stimulates the expression of AHR repressor (AHRR), 
which provides a negative feedback regulation by competing with 
ARNT for binding to AHR (Mimura et al., 1999). 

Testicular germ cell cancer derives from carcinoma in situ (CIS) 
cells, also known as intratubular germ-cell neoplasia unclassified 
(Skakkebaek etal, 1982; Oosterhuis and Looijenga, 2005), that 
are believed to arise from primordial germ cells that are blocked in 
differentiation. The initial malignant transformation of germ cells 
is thought to be caused by hormonal disturbances in the microen- 
vironment of differentiating germ cells (Skakkebaek etal, 1987; 
Sharpe and Skakkebaek, 1993; Rajpert-De Meyts etal, 1998). An 
association of TGCC with maternal estrogen and androgen lev- 
els in early pregnancy has recently been reported, supporting this 
hypothesis (Holl et al., 2009). In addition to altering transcription 
of steroidogenic genes, such as cytochrome P450 enzymes, the 
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AHR signaling pathway has been shown to interface with estrogen 
and androgen receptor signaling (Jana et al., 1999; Kizu et al., 2003; 
Safe and Wormke, 2003; Boverhof et al., 2006; Beischlag et al., 2008; 
Kollara and Brown, 2010). Additionally, AHR plays a significant 
role in tumor promotion and progression by deregulation of cell- 
cell contact and inhibition of apoptosis (Mulero-Navarro etal, 
2005; Diry et al., 2006; Carvajal-Gonzalez et al, 2009; Chopra et al., 
2009; Dietrich and Kaina, 2010). 

The risk of environment-related carcinogenesis depends not 
only on the magnitude of the exposure, but also on individual 
susceptibility to the pollutant. Therefore, we aimed to analyze the 
association between polymorphic variants of genes encoding AHR 
and AHRR and the risk as well as histological and clinical features 
ofTGCC. 

MATERIALS AND METHODS 

The cohort consisted of TGCC patients from Sweden and Den- 
mark, two neighboring Scandinavian countries, with populations 
sharing largely the same genetic background as both descended 
from a common ancestral population. 

SWEDISH TGCC PATIENTS 

All TGCC patients that were referred to the Department of 
Oncology, Skane University Hospital, Lund, the Department of 
Oncology, Radiumhemmet or Sodersjukhuset, Karolinska Univer- 
sity Hospital, Stockholm since March 1996 and November 1998 
until October 2006, respectively, were asked to participate in a 
study focusing on reproductive function. In total 460 patients 
were eligible for the study, of which 75 declined to participate 
and 45 were excluded due to linguistic problems, bilateral testicu- 
lar cancer, physical handicap, or moving to another region. Seven 
patients were excluded due to compromising mental conditions, 
10 were excluded due to contra-lateral testicular cancer diagnosed 
after inclusion and 3 died of progressive disease before blood sam- 
ples were obtained. Among the 320 participants in the study on 
reproductive function, no DNA was available for 39. Additionally, 
3 patients were excluded due to extra-gonadal germ cell cancer, 
leaving the total at 278. All patients gave their informed written 
consent and the study was approved by the Ethical Boards of Lund 
University and Karolinska Institute. 

DANISH TGCC PATIENTS 

In the period between 1999 and 2008, DNA was collected from 
approximately 300 TGCC patients on the occasion of semen bank- 
ing prior to surgery or fertility assessment after treatment. One 
hundred samples were random selected for this study; the criterion 
for selection being sufficient DNA amount. In total, 1 1 patients 
were excluded because of mistaken diagnosis ( n = 6), purely extra- 
gonadal tumor (n = 4), and in one case because the genetic SNP 
analysis failed. Among the remaining 89 patients, 5 presented with 
CIS only and were therefore not included in the analysis of distri- 
bution of seminomas (SE) versus non-seminomas (NSE). For six 
patients no information on stage was available. 

CONTROL SUBJECTS 

Control subjects were recruited in a study of reproductive func- 
tion among Swedish military conscripts aged 18-20 years in the 



period 2000-2001 (Richthoff et al., 2002). As part of the investiga- 
tion, scrotal palpation and ultrasound was performed in order to 
exclude testicular tumors or microcalcifications, which are indica- 
tive of an increased risk of CIS. Furthermore, they delivered one 
ejaculate for semen analysis as well as a blood sample for assess- 
ment of hormone levels. Among the 305 men that were eligible, 
214 men with a Swedish mother were selected as control subjects 
(Table 1). 

GEN0TYPING 

Genomic DNA was prepared from peripheral leukocytes using 
QIAamp DNA Maxi Kit (Qiagen, Germany). All samples 
were normalized to the same DNA concentration and the 
genotypes were determined using the Sequenom MassArray 
MALDI-TOF mass spectrometry. Eleven SNPs in the gene 
encoding aryl hydrocarbon receptor (AHR) and 18 SNPs in 
the gene encoding aryl hydrocarbon receptor repressor (AHRR) 
with a minor allele frequency >0.05 that were identified as 
haplotype-tagging SNPs were selected using dbSNP (available at: 
http://www.ncbi.nlm.nih.gov/SNP) and SNP assays were designed 
using MassArray Assay Design ver. 2 software (Sequenom Inc., 
USA). Primers were obtained from Metabion GmbH (Germany) 
and all reactions were run under the same conditions, except for 
the primer annealing temperature of the primary PCR. PCR reac- 
tions were performed in a total volume of 6 u,L containing 2.5 ng 
template DNA, 1.25x Taq PCR buffer (Hotstar, Qiagen), 0.15 U 
Taq polymerase (Hotstar, Qiagen), 3.5 nM MgCl 2 , 0.5 mM dNTP, 
and 100 nM of each primer. Amplifications were performed using 



Table 1 I Characteristics of controls and case patients. 


Control 


TGCC-SE 


TGCC-DK 


n (%) 


n (%) 


n (%) 


n 214 


278 


89 


Age (year) 18 ±0 


31 ±7 


31 ± 8 


Histology 


NSE 


153 (55.0) 


34 (38.2) 


SE 


125 (45.0) 


50 (56.2) 


CIS only 




5 (5.6) 


Stage 


Localized 


200 (71.9) 


65 (73.0) 


Metastasis 


78 (28.1) 


18 (20.2) 


Unknown 




6 (6.7) 


Family history ofTGCC 


Yes 0 (0) 


9 (3.2) 


7 (2.3) 


No 214(100) 


250 (90.3) 


77 (88.5) 


Unknown 


19 (6.5) 


5 (9.2) 


History of cryptorchidism 


Yes 7 (3.3) 


23 (8.3) 


13 (14.6) 


No 207 (96.7) 


247 (88.9) 


63 (70.8) 


Unknown 


8 (2.9) 


13 (14.6) 
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GeneAmp 9700 machines with dual-384 heads as follows: 95°C 
for 15 min, 45 cycles at 95° C for 20 s, 56° C, 60° C, or 64° C for 
30 s, 72° C for 60 s, and finally 72° C for 3 min. Dephosphoryla- 
tion of unincorporated dNTP was achieved using shrimp alkaline 
phosphatase. Concentrations of individual homogenous MassEX- 
TEND (hME) primers were adjusted to even out peak heights in 
the mass spectrum. The extension reactions were carried out by 
mixing the adjusted primer mix (containing approximately 1 u,M 
of each primer) with hME mix containing buffer and 50 (xM of 
each d/ddNTP mix and 1.25 U of Thermo Sequenase (Amersham 
Biosciences, Uppsala, Sweden). PCR amplification of hME reac- 
tions was performed as follows: 94° C for 2 min and 99 cycles 
at 94° C for 5 s, 52° C for 5 s, and 72° C for 5 s. The samples 
were then manually desalted by using 6 mg of Clean Resin and 
a dimple plate and subsequently transferred to a 384-well Spec- 
troCHIP using a nanodispenser. Randomly selected samples of 
each genotype were directly sequenced in order to validate the SNP 
assays. 

STATISTICAL ANALYSIS 

Single-nucleotide polymorphisms data was processed and ana- 
lyzed using the web-based SNPator data analysis suite (Morcillo- 
Suarez et al., 2008). Agreement with Hardy- Weinberg equilibrium 
was tested using a x 2 goodness-of-fit test. TGCC patients were 
divided in two groups depending on whether they were diagnosed 
with SE or NSE. Logistic regression association was used to calcu- 
late the odds ratio (OR) and 95% confident interval for developing 
TGCC, for developing either SE or NSE, and for developing dis- 
seminated disease [defined by stages II, III, and IV according to 
the Royal Marsden Hospital staging classification (Husband and 
Koh, 2004), with stage I disease defined as tumor confined to 
the testis, with no evidence of metastases] . In the latter analysis 
histological subtype was included as a covariant in the logistic 
regression. Associations were considered statistically significant at 
p < 0.05. Assuming a co-dominant additive model, we tested a 
linear trend of increasing effect in the different genotypes using 
linear-by-linear association x 2 statistics (SPSS ver. 21). Linkage 
disequilibrium (LD; i.e., non-random association of alleles) was 
assessed for each gene. The LD between identified SNPs was deter- 
mined by pairwise comparisons of correlation coefficient between 
SNPs (r) and the likelihood that recombination has occurred 
between SNPs (£>'). Since we expected that the selected SNPs 
have a small contribution to TGCC and we regard our study to 
be exploratory rather than confirmatory, we have chosen to avoid 
correction for multiple testing. This may increase the risk of type 
1 errors, but it prevents type 2 errors. 

RESULTS 

No significant deviations from Hardy- Weinberg equilibrium were 
detected for the SNPs included in the analysis. 

We did not find a statistically significant association between 
the studied polymorphisms and the risk of TGCC. 

For four SNPs in AHRR a significant allele association with 
the occurrence of disseminated disease was observed (Table 2). 
Patients with metastatic disease had significantly lower frequen- 
cies of the minor rs2466287 G, rs2672725 G, rs6879758 C, and 
rs6896163 G alleles compared to patients with localized disease 



(5 vs 10%, 6 vs 11%, 2 vs 5%, and 2 vs 6%, respectively), which 
were associated with 57, 51, 73, and 66%, reduced per allele OR 
for developing metastatic TGCC, respectively. Heterozygous car- 
riers of AHRR variants rs6879758 and rs6896163 had a 74 and 
67% reduced risk of developing metastatic TGCC, respectively. 
Due to very low frequencies, or even absence, of homozygous car- 
riers of the AHRR variants the OR could not be calculated for 
these groups of patients. However, trend analyses demonstrated 
significant associations for all four AHRR variants. 

Polymorphisms in AHRR did not show a significant association 
with risk of TGCC or histological subtype, which was also the 
case for SNPs in AHR. Histological subtype and the occurrence 
of metastasic disease were not associated with SNPs in AHR. LD 
analysis demonstrated a high correlation between rs2466287 and 
rs2672725 as well as between rs6879758 and rs6896163 (Table 3). 

When analyzed in the two populations individually, similar 
statistically significant ORs for the four SNPs were observed in the 
Swedish population as compared to the combined population. The 
associations in the Danish populations showed similar tendencies, 
although these did not reach statistical significance. 

DISCUSSION 

In this study we have analyzed associations between SNPs in genes 
encoding AHR and AHRR and the risk of developing TGCC, his- 
tological subtype, and the occurrence of metastasis. Whereas no 
associations were found with SNPs in AHR, four variants in AHRR 
associated significantly with the occurrence of metastatic disease. 

Both AHR and its negative regulator AHRR are ligand-activated 
transcription factors that belong to the family of basic helix-loop- 
helix Per-Arnt-Sim (bHLH-PAS) proteins. They are ubiquitously 
expressed in almost all human tissues, with levels particularly 
high in the testis (Burbach etal, 1992; Tsuchiya etal, 2003). 
2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is the most well- 
characterized exogenous ligand of AHR and it is known for its 
anti-estrogenic properties such as inhibition of estradiol-induced 
uterine weight increase and decreased levels of estrogen and 
progesterone receptors in uterus of rat and mice (Safe, 1995). 

Several epidemiological studies have shown that TCDD pro- 
motes cancer, and experiments in AHR-deficient mice have 
shown that AHR is essential for the tumor-promoting effects of 
TCDD (Andersson etal., 2002; McMillan and Bradfield, 2007). 
Conversely, transgenic mice with a constitutively active AHR 
spontaneously develop tumors, their weights of testis and ven- 
tral prostate are decreased and the epididymal sperm reserve is 
reduced (Brunnberg etal., 2011). Furthermore, AHRR, the nega- 
tive regulator of AHR signaling, functions as a tumor suppressor 
in multiple human tumors (Zudaire etal., 2008; Li etal, 2012). 
Several studies have shown that the AHR activation also inter- 
feres with reproductive health. An increased incidence of TGCC 
has been observed in workers exposed to TCDD (Kogevinas et al, 
1997) and dioxins have been reported to reduce sperm number, 
and to decrease accessory sex gland weight and anogenital distance 
in rats (Mably etal, 1992; Gray etal, 1995; Sommer etal., 1996; 
Faqi etal, 1998; Ohsako etal., 2002). 

While SNPs in either AHR or AHRR did not associate with the 
risk of TGCC, it is interesting that an association was observed 
between AHRR and dissemination of TGCC, since AHR signaling 
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Table 3 | Linkage disequilibrium (D ) and correlation coefficients (r) 
between SNPs in AHRR that associate with metastatic TGCC in the 
total population. 



rs2466287 rs2672725 rs6879758 





D' 


r 


D' 


r 


D' r 


rs2466287 
rs2672725 
rs6879758 
rs6896163 


0.978 
0.239 
0.195 


0.909 
0.172 
0.142 


0.230 
0.186 


0.155 
0.132 


0.999 0.948 



is indeed known to contribute to the control of cell adhesion 
and migration (Carvajal-Gonzalez et al., 2009; Dietrich and Kaina, 
20 1 0) . A potential role of AHR signaling in dissemination of TGCC 
is further supported by Koliopanos etal. (2002) who reported 
that AHR agonists decreased anchorage-independent pancreatic 
cancer cell growth. 

Since testicular cancer rarely occurs before the onset of puberty, 
the malignant transformation of CIS cells may be associated 
with activation of the hypothalamic-pituitary-gonadal axis. An 
interaction between AHR signaling and steroid or gonadotropic 
hormones could therefore be an alternative mechanism underlying 
the association between AHRR polymorphisms and metastatic 
TGCC. 

A strength of our study is the access to two independent TGCC 
populations. Although not all findings reached statistical signif- 
icance in the Danish cohort when the countries were analyzed 
separately, most probably due to the limited number of cases 
in this group, the same trends with similar relative risks were 



found, which strengthens the relevance of our findings. We cannot 
exclude that a small fraction of the TGCC patients in our cohort 
were from non-Caucasian background. However, none of the allele 
frequencies of the analysed SNPs deviated from Hardy- Weinberg 
equilibrium, which indicates that either the non-Caucasian allele 
frequencies did not differ from those in the Caucasian popula- 
tion, or that non-Caucasians did not significantly contribute to the 
total cohort. The SNPs identified in this study are intronic poly- 
morphisms, thus excluding structural changes in the proteins that 
these genes encode. Yet, polymorphisms in non-coding regions 
can still play an important role in the regulation of gene expres- 
sion, e.g. , by affecting mRNA stability, regulation by intronic micro 
RNAs, or potential splice sites. The identified SNPs could also be 
in LD with other genetic polymorphisms that are causally related 
to TGCC. 

In summary, we have identified polymorphic variants in the 
gene encoding AHRR that associated with risk of disseminated 
TGCC. This supports our hypothesis that the AHR signaling path- 
way may affect the progression of TGCC, possibly by interfering 
with reproductive hormone actions. 
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